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ABSTRACT 
(N . 

^D ' We study a general elliptical potential of the form ^p{x^ + y'^/q^) {0 < q < 

f^ ' 1) plus an additional shear (with an arbitrary direction) as models for the 

Qs , observed quadruple lenses. It is shown that a minimum additional shear is 



needed even just to reproduce the observed positions alone. We also obtain 

the dependence of the axial ratio, q, on the orientation of the major axis of 

O 

Jh I potential. A general relation also exists between the shear, the position angle 

jrt . and axial ratio of the lensing galaxy. The relation shows a generic degeneracy 

in modelling quadruple lenses. In particular, it shows that only the ratio of the 
• 1—1 . 

^^ ' ellipticity, e = {l—q-^)/{l+q'^), to the magnitude of shear, 7 can be determined. 

Vh ' 

C^ ' All these results are valid regardless of the radial profile of the potential. Our 

formalism applies when the galaxy position is observed, which is the case 

for seven of the eight known quadruple lenses. Application to these seven 

cases reveals two quadruple lenses CLASS 1608+656 and HST 12531-2914, 

requiring highly significant shear with magnitude w 0.2. For HST 12531-2914, 

there must be a misalignment between the major axis of light and the major 

axis of potential (mass). We conclude that detailed modelling of quadruple 

lenses can yield valuable quantitative information about the shape of lensing 

galaxies and their dark matter halos. 

Key words: galaxies: structure ~ gravitational lensing - quasars: individual 
(CLASS 1608+656, HST 12531-2914) 
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1 INTRODUCTION 

Up to now, more than twenty multiply imaged quasars have been discovered, with roughly 
half double and half quadruple lenses (Keeton & Kochanek 1996, hereafter KK96, see also 
Schneider, Ehlers & Falco 1992, Blandford & Narayan 1992 and Kochanek & Hewitt 1996 
for general reviews). These systems, in particular the quadruple lenses, provide a unique tool 
to probe the potentials of galaxies (e.g., Kochanek 1991; Wambsganss & Paczyhski 1994; 
Witt, Mao & Schechter 1995). Recently, Keeton, Kochanek & Seljak (1996) found that while 
simple models such as singular isothermal density ellipsoids provide a reasonable statistical 
model for the whole sample, no individual quadruple system is well fitted by such models. 
They found numerically that an additional shear term can drastically reduce x^ in fitting 
while changing the radial distribution of the potential helps little. As only a limited number 
of radial profiles (such as power laws) have been explored numerically, it is not clear whether 
the bad fits can still be due to our incomplete knowledge of the radial profile of galaxies. In 
this paper, we study a general class of elliptical potentials of the form ipi^x^ + y'^/q'^), where q 
denotes the axial ratio of the elliptical potential. We show that a minimum additional shear 
is required even if one is trying to fit only the observed positions. We also obtain analytical 
formulae for the axial ratio and orientations of the potential and shear. This implies some 
generic degeneracy in modelling of quadruple lenses. Our formalism applies regardless of the 
functional form oi ip, as long as the lensing galaxy position is observed. 

The results presented here (§2) complement the study by Keeton et al. (1996) and provide 
an analytical understanding of their results. Our analytical formalism makes it possible to 
check quickly whether any elliptical potential or density distribution without shear can work 
at all without computing x^- We apply the formalism to seven of the eight known quadruple 
lenses, including three of the four quadruple lenses studied by Keeton et al. (1996). Our 
results are consistent with theirs for these three. For two of the other four cases, HST 
12531-2914 and CLASS 1608+656, we found that the minimum shear required is ^ 0.2. For 
HST 12531-2914, the major axis of potential must be misaligned with that of the light. The 
origin of the large additional shear and its implications are discussed in the last section. 



* E-mail: hwitt@aip.de 

T E-mail: smao@mpa-garching.mpg.de 



© 1997 RAS, MNRAS 000, §-0 



2 ELLIPTICAL POTENTIAL PLUS SHEAR 

Both the elhptical density and potential distributions are widely used in the literature to 
model gravitational lenses (e.g., Blandford & Kochanek 1987, Kochanek & Blandford 1987, 
Kormann, Schneider, & Bartelmann 1994a,b). These two distributions resemble each other 
when the ellipticity is small (Kassiola & Kovner 1993), which is generally the case for the 
known quadruple lenses. We will use the elliptical potential model due to its analytical 
simplicity. Our results should apply to the elliptical density distribution almost equally well. 
To be quadruply lensed by an elliptical potential, the source must be located closely behind 
the centre of the lensing galaxy. For a pure elliptical potential, the possible locations of the 
images and the lensing galaxy are highly restricted (Witt 1996, hereafter W96). In reality, 
the image positions and flux ratios will depart from those predicted by a pure elliptical 
potential. For example, large scale structure and/or other galaxies along the line of sight 
can distort the image configuration. In addition, any departure of the galactic potential from 
the idealized elliptical form can produce deviations as well. As the pure elliptical potential 
can reproduce the overall observed image configuration quite well, a reasonable approach 
for further refinement is to model all the other perturbations as an additional shear term in 
the lens equation (Kovner 1987). This is the approach we will adopt here, as in Keeton et 
al. (1996). 

We therefore assume the potential can be modelled as a two-dimensional elliptical poten- 
tial plus an additional shear in an arbitrary direction. The elliptical potential is by definition 
given by ip{re), where Ve = x"^ + y'^/q'^, and g (0 < g < 1) is the axial ratio of the potential. 
The centre of the galaxy is always located at the origin. For clarity, we first assume that 
the X-axis coincides with the major axis of the lensing potential. The results derived in this 
special coordinate system (which we term as the major axis frame) are then generalized to 
the case with an arbitrary major axis orientation afterwards. Throughout the paper, all the 
quantities measured in a general coordinate system will have a prime superscript to avoid 
confusion with those measured in the major axis frame. 

The (projected) surface mass distribution is given by Aip = 2K{x,y), where K{x,y) = 
S(x, |/)/Scrit is expressed in units of the critical surface mass density Scrit which depends on 
the distances to the defiector and the source (cf. Schneider et al. 1992). The two-dimensional 
defiection angle is then simply given by the derivatives of the potential, cx.='ViIj plus the two 
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terms related to the shear. The lens equation can be written as 

^ = a; + 7ia; + 722/ ^^ = x + 7ix + 722/ 2x, (1) 

ax are 

d^{x'^ + y"^ / q^) di){re)2y 

V = y + l2X--iiy ^ = y + l2X-'^iy ^ 5-, (2) 

ay OTe q 

where (^, 77) denotes the (unknown) source position and the magnitude of the shear is 



given by 7 = \Jl'l + l2- The shear can also be written in a "vector" form (71,72) = 
(7 cos 26*^, 7 sin 26*^) ( < 0^ < tt). Notice that the factor of 2 before 0^ arises because 
the shear is not really a vector but a tensor. When the shear is acting on-axis, we have 
9^ = 0,72 = 0. The shear is maximum off-axis when 71 = 0, i.e., when 9^ = 45°, or 135°. 

For quadruple lenses, the positions of the four images obey eqs. (|I]) and (0), therefore 
for each of the four images we can eliminate the factor dijj{re)/dre to obtain the following 
equation: 

= q — tor 2 = 1,.. .,4. (3j 

V-Vi- l2Xi + 'JWi yi 

In the next three subsections, we will use eq. (P) as basis to derive analytical results in this 
paper. 

2.1 Lower Limit On the Additional Shear 

Using the four equations as in eq.(H), we can eliminate first q and then ^ and 77, which leads 
us to the following equation: 

7iai + 72^2 + as = 0, (4) 

where the coefficients ai, a2 and 03 depend on the four relative image positions, and are 
given by 

ai = (xi + yf)f 234 - (xl + yl)f 341 + (xl + 2/3) /412 - (xl + yl)f 123, (5) 

0-2 = -yfh234 + ylhii - ylhii2 + ylhi23, (6) 

as = {xl - yl)f234 - {xj - yl)f34i + (3:3 - yl)f 412 - {xl - yl)fl23, (7) 
with the functions fijk and hijk defined as 

fijk = Xiyi[xjyk - Xkyj] + Xjyj[xkyi - xiyt] + Xkyk[xiyj - Xjyi] (8) 

hijk = xl [xjyk - XkVj] + x] [xkVi - Xiyk] + xl [xiyj - Xjyi] . (9) 

Note that fijk-, hijk, ^i; 02; and a^ are all odd under the permutation of any two indices. 
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When we derived eq. (^ we assumed that the major axis of the lensing potential is along 
the X-axis. In practice, it is usually more difficult to measure the orientation of the galaxy 
than to measure its position and those of the images. In any case, what we measure is the 
position angle of the light distribution, not that of the potential (mass) which enters the lens 
equation. Hence, we need to generalize eq. (^ to the case when the major axis of the galaxy 
potential is unknown, i.e., we need to investigate what happens to the coefficients oi, 02 and 
03 when the major axis of the lensing galaxy is rotated by an angle ^g(— 7r/2 < 6g ^ ^/2)- 
The positions of the images in the new coordinate system, {Xi,y'j), are related to those 
measured in the major axis coordinate system, {xi,yi), by 

Xi = x'iCos9G + yiSm9G, (10) 

yi = -x'i sin ec + y'i cos 9g, (11) 

for i = 1, ..,4. 

In the appendix we show that 03 is rotationally invariant, i.e., 

03 = 4. (12) 

Since 6^ = 6' — 6g, the shear tensor transforms like 

71 = 7 cos(20^) = j[ cos{29g) + 72 sin(2^G), (13) 

72 =7sin(2^^) = -7;sin(2^G)+72Cos(2eG). (14) 

To satisfy the invariance of 03, (01,02) must transform like a tensor as well. It is easy to 
verify that (01,02) transforms as follows 

oi = OiCos(26'g) + 02sin(26'G), (15) 

02 = -o'isin(2^G) + a2Cos(2^G)- (16) 

Substituting eqs. (|T2l), (|T5|), and ([T6|) into eq. (||), one obtains 

7io'i + 72«2 + 4 = 0. (17) 

The equation has the same form as eq. (^, but now all the quantities are evaluated in the 
general coordinate system. Replacing the two shear components by 7^ = 7cos(26'^) and 
72 = 7 sm{29'^) yields 

7K cos(2e;) + o^ sin(2^;)] + o^, = 0. (18) 

The above equation can be rewritten as 



7Vo'i' + o'2' sin(2^; + (/?') + a'^ = 0, (19) 



7 
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where ip' {—n < ip' < n) is the polar angle of the vector (a'j, a[). Eq. (|1^) implies that, to fit 
the observed positions, a minimum shear is required: 

7mm = , f^ ^ = . /, ., , when 6' = 6' -^^ = -%: - jsignlag) + /cvr, (20) 

Vaf + a'2 Vai' + aa^ 2 4 

where throughout the paper k is an integer that make the angle at the left hand side of the 

equation {6'^ here) fall into the right range. 

We now make some general remarks about eq. (ppj). First, the minimum shear is required 

no matter what the radial profile is as long as the iso-potential contours are ellipses. This 

easily explains why changing radial profiles, such as adding a core radius or changing the 

slope of a power law radial profile, will not improve the fitting much (Kochanek 1991; 

Wambsganss & Paczyhski 1994; Keeton et al. 1996). Second, as the coefficients 01,02 and 

03 involve differences of permuting terms, the required (minimum) shear is likely to be very 

sensitive to the accuracy of positions. High quality relative astrometry of galaxy and image 

positions are thus much desirable. 

2.2 Axial Ratio of Lensing Galaxy 

In this subsection we discuss whether the axial ratio q can be restricted. To do this, we again 
start with eq.(^ and eliminate successively ^, 77 using three image coordinates i = 1,2,3, 
after which we obtain 

(l + 7i-g^(l-7i))/i23 = 'y2[q'^hi23-ylix2y3-X3y2)-ylix3yi-xiy3)-yi{xiy2-X2yi)].i2l) 

where /132 and hi32 are defined as in eqs. (||) and (P). If 72 7^ 0, then we can use the fourth 

image position to eliminate 71 or 72. By eliminating one shear component the other shear 

component factorizes out of the equation simultaneously. Therefore we obtain an equation 

which depends only on the relative image positions and the axial ratio: 

2 ^ Z/?/234 - yi/341 + 1/3/412 - ^4/123 ^ Q1-Q3 ^221 

3^1/234 ~ 2:2/341 + 2:3/412 — 3^4/123 (^i + % 
where Oi and 03 are given by eqs. (|^) and (|^. We can easily generalize eq. ( p2D to an arbitrary 

major axis orientation by using eqs. (pISJ) and (0): 
2 _ o-'i cos(26'g') + a'2 sin(26'G) - g^ _ sin(v9^ + 26'g) - sign(4)7min 
a'l cos(26'g) + a'2 sin(26'G) + 4 sin(v9' + 29g) + sign(a^)7i„in ' 
with ip' as defined below eq. (0). Since g^ must be positive, eq. (|23| ) provides a general (but 

weak) constraint on the axis orientation (cf. Fig. 3). Clearly q achieves a maximum, 

gmax = -- — ^^ , when Og = O^ra^ = -— + -sign(a3) + k-K. (24) 

V 1 + 7min / 2 4 
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As shown in W96, it is impossible to determine q for a pure elliptical potential with a 
shear acting along the axis, i.e., when 72 = 0; it therefore seems rather peculiar that for the 
more general shear case it is actually possible to do so. The reason for this peculiarity is 
rooted in the properties of an pure elliptical potential. Due to its highly symmetric and self- 
similar shape of the iso-potential contours, the positions of the images and galaxy are highly 
restricted. For example, the location of the lensing galaxy is restricted to a hyperbola-like 
curve and the image positions in such a potential must fulfil the identities /123 = /134 = 
/234 = /i24 = (W96). This can be seen from eq. (|2T|) by setting 72 = 0. Then we must 
have [1 + 7i — g^(l — 7i)]/i32 = 0, which requires either g^ = (1 + 7i)/(l — 7i),7i < or 
/i32 = 0. If the former condition is satisfied, using eqs. (1) and (2), it is easy to show that 
all the images must lie on a straight line. This image configuration clearly does not resemble 
the observed quadruple lenses. Therefore we must have /123 = 0. As a result q can no longer 
be determined. Physically we can understand it as follows: A off-axis shear on top of the 
elliptical potential breaks down the high symmetry required for the image positions {fijkS 
are no longer required to be zero), which in turn allows us to determine the axial ratio. 

As we have shown, when the shear is acting on-axis, eq. (|23| ) cannot be applied. Nu- 
merically this implies that when the shear component (72) is small, eq. ( PHj ) is likely to be 
unstable due to the errors in the image and galaxy positions. Therefore before applying eq. 
(^) it is necessary to check first whether the observed system requires an off-axis shear, 
using the test introduced by W96 (see §2.4). If a significant off-axis shear is indicated, the 
(minimum) magnitude of the additional shear can be estimated using eq. (^Ol) . If the required 
shear is large, then eq. (^31) can be applied to restrict the axial ratio and the orientation of 
the major axis of potential (see below). 

2.3 Orientations of Lensing Galaxy and Shear 

Combining eqs. ([191) and (^), we arrive at a new relation: 

/„/! ,^ sin26'^ 1 — 0^ _^^ 

tan2^G + ¥'' = ^^-V' ^ = T— ^' 25 

cos 20^ + 6/7 1 + q'^ 

where e as defined is the ellipticity of the iso-potential contours. This equation bears some 
similarity to eq. (22) found by Keeton et al. (1996). 

Eq. ( P^D shows two important points. First, since only the ratio e and 7 enters eq. (P3[), 
there is a degeneracy between these two parameters. This can be understood as follows: an 
increase in ellipticity (0 < e < 1) (a decrease of axial ratio of q) stretches the images more 
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along the x axis, while an increase in the shear stretches the images more along the y axis, 
the balance between these two competing effects introduces the degeneracy. Eq. ( pS]) also 
shows that even if the ellipticity and the galaxy orientation are known (e.g., if we use those 
for the light distribution), the shear tensor still cannot be determined uniquely (see Figs. 8 
and 9 in Keeton et al. 1996). 

If the potential is purely elliptical, i.e., 6^ = (72 = 0), then from eq. (|25|) , we have 

^G,purc = -y for 72 = 0. (26) 

The above expression can be shown to be identical to eq. (7) in W96 in this case. However, 
if an off-axis shear is present (0 < 72 ^ 0.3), eq. (|26|) is only approximately true. Combined 
with eq. (^41), we obtain 



71" 
l^'cmax — ^'cpurcl ~ J- (27) 

Notice that in general, the difference between the true orientation and that obtained by 
using a pure elliptical potential depend on not only the direction of shear but also the ratio 
of ellipticity and shear as well. 

2.4 Application To Known Quadruple Lenses 

Eight quadruple lenses are known (see KK96 for a thorough summary). Seven of these eight 
systems (except H1413+117) have known galaxy and image positions; our results can applied 
to study these systems. When multiple sets of positions are available, we generally took the 
data with the best astrometry. Detailed references are listed in the last column in Table 1. 
For each system, we use eq. (|20|) to derive the minimum shear; the results are shown in Table 
1. As one can see, the required minimum shears vary significantly from system to system. 
For example, for 2237+0305 it is consistent with zero, while for the lens CLASS 1608+656 it 
is as large as 0.25. How sensitive are these estimates to the positional errors? To address this 
question, we used Monte Carlo simulations to generate synthetic lensed systems by assuming 
all the positions are uncorrelated and their errors are Gaussian. For each lens, 10,000 Monte 
Carlo realizations are generated, and the minimum shear is calculated for each of these. The 
average and standard deviation are then computed and shown in Table 1. In some cases, 
the distributions of the inferred minimum shear are highly non-Gaussian, especially in the 
case of B1422+231. It has an almost fiat distribution of minimum shear between to ~ 0.3. 
For MG 0414+0534 the minimum shear is about 0.1 (2cr significant). However, it remains 
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Table 1. Minimum Shear Required For Known Quadruple Lenses 



Object 


Tmin 


^imagc 


f^galaxy 


Reference 


2237+0305 


0.0092 (0.0099 ± 0.0063) 


0.005 


0.005 


C91 


PG 1115+080 


0.053 (0.058 ± 0.034) 


0.005 


0.05 


K93 


MG 0414+0534 


0.12 (0.12 + 0.058) 


0.0003 


0.05 


KMH96; F96 


CLASS 1608+656 


0.25 (0.25 + 0.027) 


0.01 


0.01 


M95; S95 


B 1422+231 


0.11 (0.165 + 0.125) 


0.002 


0.05 


P92; YE94 


HST 14176+5226 


0.037 (0.040 + 0.027) 


0.03 


0.03 


R95; R96 


HST 12531-2914 


0.18 (0.18 + 0.058) 


0.03 


0.03 


R95; R96 



Note. — 7min is the minimum shear which is required to fit the positions of the observed image and lensing galaxy. The mean 
and standard deviation (listed in brackets) for each system in 7inin E're derived from 10,000 Monte Carlo simulations for the 
positions assuming the positional errors are Gaussian and uncorrelated. criini<.gc and o-gaiaxy denotes the astrometric accuracy 
(in arcseconds) of the image positions and the lensing galaxy (cf. KK96). One known quadruple lens H1413+117 is not listed 
here because no galaxy position is yet available, therefore our formalism does not apply. 

unclear whether the shear is due to the less accurate galaxy position ^. Two cases, CLASS 
1608+656 (Meyers et al. 1995, hereafter M95) and HST 12531-2914 (Ratnatunga et al. 1995, 
1996, hereafter R95, R96) require even more significant shears than MG 0414+0534. We will 
therefore concentrate on these two systems below. 

The image configurations of these two systems, HST 12531-2914 and CLASS 1608+656 
are shown in Figs. 1 and 2, respectively. For a pure elliptical potential, the lensing galaxy 
position must lie on the hyperbolic-like curves in these figures. Both systems show a signifi- 
cant deviation these curves, suggesting the presence of a significant off-axis shear (cf. W96). 
This is consistent with the highly significant minimum shear derived above (cf. Table 1). 

As HST 12531-2914 and CLASS 1608+656 clearly need some large off-axis additional 
shear, eq. (P^D can therefore be applied to both systems. In Fig. 3, we plot the dependence 
of the axial ratio q on the position angle 6p for these two systems. Note the position angle 
is conventionally measured from north (positive y axis) through east (minus x axis). If we 
limit the range of 9p to be from to tt, then it is related to 6g by 

Op = 60 + 7r/2. (28) 

If we (artificially) fit a pure elliptical potential model to both systems, we obtain position 
angles Op = 18.6°, 9p = 68.5° for HST 12531-2914 and CLASS 1608+656, respectively (cf. 

T The EUithorpe (1995) image position (sec KK96) would imply a considerable minimum shear (fa 0.16) for MG 0414+0534. 
However, his relative position of image C is not compatible with the recent observations of Falco et al. (1996) and Katz et al. 
(1996). 
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Figure 1. Image configuration of HST 12531—2914. North is up and east is to the left. The errorbar in the galaxy position is 
about 0.03 arcsecond (cf. Table 1). The solid line indicates th possible location of the lensing galaxy when we assume a pure 
elliptical potential. For this model, the galaxy and image positions lie on the same curves. The offset of the galaxy position from 
the predicted curves indicates the presence of a large shear. The expected position angle, dp, measured from north through east, 
for a pure elliptical potential is indicated with a dashed line (cf. W96). The observed position angle for the light distribution 
is 22° .6 ± 0.5 from the HST image in the F606W filter (R95; R96). 

Fig. 3). For HST 12531-2914, this estimate is in rough agreement with the position angle of 
the hght, Op = 22°.6 ± 0.5 (R95). (For CLASS 1608+656, the orientation of the galaxy is 
unavailable.) There are a few interesting things that can be seen from Fig. 3. First for HST 
12531-2914, the positional angle of light (solid line) falls in the unphysical region where q^ 
is negative. This means that there must be a misalignment between the major axis of the 
potential and that of the light. This conclusion is valid regardless of the radial profile of the 
potential. From Fig. 3, if the axial ratio of the potential is identical to that of light (0.73), 
then 6p ^ 39.5°, implying a misalignment of about 17°. Second, the q vs. 6g curves are fairly 
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Figure 2. Image configuration of CLASS 1608+656. Nortti is up and east is to tlie left. Ttie errorbar in the galaxy position is 
about 0.01 arcsecond (cf. Table 1). The solid line indicates the possible location of the galaxy when we assume a pure elliptical 
potential. For this model, the galaxy and image positions lie on the same curves. The offset of the galaxy position from the 
predicted curves indicates the presence of a large shear. The expected position angle, 6p, for a pure elliptical potential is 
indicated as a dashed line (cf. W96). The inferred position angle from modelling is 67° (E of N) (M95). 



flat around the peak regions. This can be seen by by expanding eq. (|2J) in Taylor series 
around 6'G,max5 which leads to 

?' = dax - TT^T^'^C ^G = ^G,max + 5^G. (29) 

1,-L -r /minj 

This means that q^ varies quadratically at the peak region. Third, although we showed that 
HST 12531-2914 and CLASS 1608+656 cannot be exactly fitted by a pure elliptical potential, 
nevertheless the models published so far did use the pure elliptical density distributions (R95, 
R96; M95). In these modelling, one finds the best fit parameters by minimizing a x^ measure. 
Obviously the resulting x^ per degree of freedom will be very large (see x^ in Keeton et al. 
1996 for other systems). As one typically searches for the best fit axial ratio starting from an 
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initial guess of one, as shown in Fig. 3 the q vs. Oq curve is nearly flat around the maxinium 
(implying a large phase space in the multiple dimensional parameter space), we therefore 
expect the resulting axial ratio to be close to the maximum axial ratio . Indeed, the axial 
ratios for the density distribution, gp, are found to be 0.37 for HST 12531-2914 (R95; R96) 
and 0.28 for CLASS 1608+656 (M95). Since g ^ 2/3 + qj?, (Binney & Treimaine 1987), 
we have q = 0.83 and 0.79. These are very close to the maximum axial ratios from Fig. 
3, gmax = 0.83 for HST 12531-2914 and g^^x = 0.78 for CLASS 1608+656, achieved at 
6p = 65.3° and 114.2°, respectively. Fourth, the position angle at the maximum axial ratio 
is approximately 45° away from the orientation inferred by modelling the lens galaxy as a 



pure elliptical potential for both systems, just as given by eq. (|27|) 



3 DISCUSSION 

We have studied a general class of models with an elliptical potential plus an additional 
shear. It was shown that to fit the image and galaxy positions in quadruple lenses, the 
additional shear has to exceed some minimum value. In addition, an analytical expression 
for the axial ratio is derived. We also showed that there is a complex relation between the 
orientations of the shear and potential, the ellipticity and the magnitude of shear. As only the 
ratio of the ellipticity and the magnitude of shear enters the relation, these two parameters 
are linearly degenerate. We emphasize our results are valid no matter what the radial profile 
is for the elliptical potentials. Applying the analytical results to seven of the eight known 
quadruple lenses, we found that two (MG 0414+0534, B 1422+231) are consistent with 
the presence of additional shears of the order of 0.1, while HST 12531-2914 and CLASS 
1608+656 require highly significant shears of ^ 0.2. For HST 12531-2914, we conclude that 
the major axes of potential and light must be misaligned, regardless of the detailed potential 
shape. We caution that both systems are somewhat "special": HST 12531-2914 has very 
small separations between the images. The images seem to be not perfectly aligned with 
axes of the lensing galaxy (cf. the frames in R96). In contrast the lensed source in CLASS 
1608+656 is a radio galaxy and not a quasar. It is conceivable that the small separations or 
the extended source size make the position measurements more difficult and the errors on 
their positions are under-estimated. We therefore artificially enlarged their errors by a factor 
of 2, and recomputed their statistical significance. CLASS 1608+656 remains 5a significant 
while the significance for HST 12531-2914 has dropped to 1.7cr. The exceptional nature of 
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Figure 3. Dependence of the axial ratio q on the position angle 9p for CLASS 1608+656 (bottom panels) and HST 12531-2914 
(top panels). The position angle of the light distribution for HST 12531-2914 in the F606W filter is shown with the solid line. 
The dashed line shows the position angle of the major axis estimated by enforcing a pure elliptical potential while the dotted 
line shows the position angle for the maximum q (see eq. [M). gmax = 0.78 for CLASS 1608+656 and gmax = 0.83 for HST 
12531-2914. 

CLASS 1608+656 can also be seen from the ratio of the distances between the images C and 
D to A and B, AOcd/^Oab ~ 1-0 (for HST 12531-2914 Mcd/^Oab ^0.73). For a pure 
elliptical isothermal sphere, we would expect /S.6cd/ ^Oab ~ <i- CLASS 1608+656 therefore 
signals a substantial deviation even from the maximum possible value (g = 0.78), again 
indicating the presence of an additional shear on top of an elliptical potential. 

We emphasize here that in Table 1 only the lower bound of the additional shear was 
derived. Computations with some test potentials indicates that 7min ~ 72 and gmax ~ 1 
when 7i ^ 0, i.e., when 9^ ~ 45°, or 135°. In these cases the additional shears derived are 
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highly significant. However, if 71 ^ 72 the minimum shear derived using eq. ( pO] ) is usually 
smaller than the actual shear 7 and not statistically significant. Since for CLASS 1608+656 
and HST 12531-2914, the additional shears are highly significant, their directions of shear 
are likely to be close to 45° and the shear is close to the minimum values derived here. For 
B 1422+231 and PG 1115+080, Keeton et al. (1996) gives the additional shears of 0.20 and 
0.09, roughly a factor of 2 larger than the predicted minimum shears, suggesting there is a 
significant part of shear acting along the axes. The real shears in other systems can also be 
easily a factor 2 or 3 larger than the predicted values. 

These additional shears required are large compared with the external shears produced 
by large scale structure (Bar-Kana 1996, Schneider 1997), or galaxies and clusters along 
the line of sight (Keeton et al. 1996), which are usually of the order of a few percent. This 
suggests that the origin of the shear is not "external" , but introduced by the lensing galaxy 
internally (Keeton et al. 1996). If the shear is truly intrinsic, it will be interesting to see 
whether there is correlation of the additional shear with physical parameters of the lenses. 
For example, the potential of galaxies will be less relaxed and more irregular in a hierarchical 
formation scenario as the redshift increases. Therefore when a sufficient number of lenses 
are known, one should find some correlation between the required additional shear and the 
redshift of the lensing galaxy. 

It is worthwhile to reflect why we need such large additional shears to model some 
of these observed systems. For the elliptical potential studied here (and similarly for the 
elliptical density distribution), the direction of the deflection angle is independent of the 
radial coordinates (cf. eqs. [1] and [2]). This strongly restricts the allowed image and galaxy 
positions (W96). Violation of these restrictions directly translates to an additional off-axis 
shear required in the model. The large values of the inferred additional shears illustrates 
that our modelling of the lens potential is too simplistic. The misalignment of the mass 
and luminous part of the lensing galaxy in HST 12531-2914 shows an example of possible 
complexities. Other possibilities such as the triaxiality of the dark halo clearly exist (see 
Keeton et al. 1996). All these complications make the isopotential contours more complex and 
possibly twisted. With the added complexities, one can presumably fit the observed positions 
and flux ratios better as more degrees of freedom become available. Such complexities may 
prove to be a nuisance in lens applications such as determining the Hubble constant. On the 
other hand, this means detailed modelling of quadruple lenses may also yield information 
about the shapes of lensing potentials and dark halos. For example, the newly discovered 
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quintuple lens 0024+1654 (cf. Colley, Tyson & Turner 1996) would be an excellent example 
to apply our formalism as the fifth image provides additional constraints. This exceptional 
case shows a lensed high-redshift galaxy with large image separations. Since the source is also 
extended, the system contains much more information that can be used to probe the potential 
of the lensing cluster. With more and more quadruple (or quintuple) systems discovered, we 
are optimistic that gravitational lenses will become an increasingly discriminating tool to 
study the (dark) matter distribution in galaxies. 

This work was supported by a postdoctoral grant of the Deutsche Forschungsgemein- 
schaft (DFG) under Gz. Mu 1020/3-2 (HJW) and by the "Sonderforschungsbereich 375-95 
fiir Astro- Teilchenphysik" der Deutschen Forschungsgemeinschaft (SM). We are very grateful 
to Peter Schneider for his constructive comments on the paper. 

APPENDIX A: ROTATIONAL INVARIANCE OF A, 



In this appendix, we show that a^ as defined in eq. (|^ is rotationally invariant. It is easy to 
verify that a^ can be written as the determinant of a 4x4 matrix: 
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From linear algebra, this determinant can be expressed as 

.2\ / 



03 = ^ e{ij,k,l) det 

i,j,k,l 



X det 




(A2) 



2 22 
xf - yf Xj - Vj 

where {i,j,k,iys are permutations of the four indices (1,2,3,4) satisfying i < j and k < I, 

and e{i,j, k, /)'s are either +1 or —1 but are of no importance here, since we will show each 

term in the sum is rotationally invariant. To prove this, let us express the image positions 

in the polar coordinates, i.e., {xi,yi) 

submatrices are then 
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rj cos^j,rj sin^j). The determinants of the two 2x2 
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1 »j 
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Both terms are obviously invariant under rotation. It follows that 03 is rotationally invariant 
as well. 
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